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1 Introduction 



The question of space-time noncommutativity has a long-standing story. It was put forward by Snyder 
[1], Heisenberg, Pauli [2] and Yang [3] as a means to regularize quantum field theories. However, the 
development of rcnormalization techniques and certain undesirable features of noncommutative theories 
such as the breakdown of Lorcntz invariancc have hindered further research in this direction. More recently, 
several important developments in various approaches to the quantization of gravity have revived the 
interest in the concept of noncommutative space-time. See for instance [4] in the context of 3d gravity. In 
the realm of string theory, the discovery that the low energy effective theory of a D-brane in the background 
of a Neveu-Schwarz-Ncvcu-Schwarz B field lives on a space with spatial noncommutativity [5], [6], [7] has 
triggered an enormous amount of research in this field. From another perspective, a simple heuristic 
argument, based on Heisenberg's uncertainty principle, the equivalence principle and the Schwarzschild 
metric, shows that the Planck length seems to be a lower bound on the precision of a measurement of 
position [8]. This reenforces the point of view that a new (noncommutative) geometry of space-time will 
emerge at a fundamental level. The simplest way to implement these ideas in quantum field theories is by 
adding to the phase-space noncommutativity of quantum mechanics a new space-time noncommutativity. 
Conncs' noncommutative geometry [9], [10] provides a suitable mathematical framework for the formulation 
of quantum field theories in noncommutative space-time. These theories have been intensively studied in 
the literature (for a review see [11], [12]). In most of the approaches [10], [13] the space-time coordinates 
x'^ do not commute, either in a canonical way 

[a;'^,^'^] =^0'^^ (1) 

or in a Lie-algebraic way 

[x^x^=^C^V, (2) 

where 6^^^ and C^^ arc constants. In general, one assumes that 6^^ = (or (7^* = 0) (i.e. that time is an 
ordinary commutative parameter) to avoid problems with the lack of unitarity.^ Equations (1,2) can be 
regarded as a particular case of the more general deformation [a^^jX*"] = if^^'^{x) [15, 16], where f^^{x) are 
real functions of x such that f^'^{x) = —f'^^{x)'^. This occurs when the laboratory frame has a space-time 
dependent motion. As proposed in [18] it is likely that f^'^ has a fixed value in the Cosmic Microwave 
Background Radiation frame, which may be considered as approximately fixed in the celestial sphere. For 
this and technical reasons, many authors have chosen the simpler situation (1) to explore some of the 
qualitative (and possible quantitative) effects of the noncommutative extensions. Indeed, the constant 
matrix 6^'^ in (1) defines a bilinear skew-symmetric form 0,g{a, b) = ^ a^O^^hi,. If the dimensionality is 
even and this form is non-degenerate (this will be the case in the present work), then by a linear version of 
Darboux's Theorem [19] there exists a linear transformation which takes VLq into the standard symplectic 
form of classical mechanics (see below). This simplifies dramatically the approach to such systems. For 
these reasons, in this work we shall concentrate on the canonical case (1), and leave other extensions (such 
as (2)) for future a future work. 

Adopting this point of view many authors have addressed the problem of testing the existence of space- 
time noncommutativity in nature by resorting to a quantum mechanical approximation to the complete 
quantum field theory. This theory of quantum mechanics with noncommutativity of the canonical type (1) 
is conventionally known as noncommutative quantum mechanics [20]- [31]. This sort of theory also appears 
in the context of quantum constrained systems [32] . 

Ho and Kao [26] derived quantum mechanics from noncommutative quantum field theory in the nonrcl- 
ativistic limit and showed that particles with opposite charges display opposite noncommutativity. Bolonek 

^This claim seems debatable. See [14] for an alternative. 

^see [17] for a recent survey of the various types of deformations and some of there consequences at the classical level. 
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and Kosinsky [27] proved that Heisenberg's uncertainty principle is somewhat more restrictive in noncom- 
mutative quantum mechanics in the sense that e.g. on the noncommutative plane there arc no quantum 
states saturating simultaneously more than one of the uncertainty relations. Most of the research thus 
far aims at finding upper bounds on the noncommutative parameters or at evaluating noncommutative 
corrections to ordinary quantum mechanical systems such as the harmonic oscillator, the Landau problem, 
the hydrogen atom, Lamb shift etc. [15]-[31]. Unfortunately, all these corrections seem to be too small to 
be detected experimentally with present technology. 

In this work we shall consider the plane with spatial noncommutativity: 



[qi, qj] = iOeij, [qi,Pj] = ihSij, \pi,Pj] = 0, i,j = 1, 2, 



(3) 



where ei2 = —£21 = 1, en = £22 = and is a real constant. We shall call this the extended Heisen- 
berg algebra. All the particles appearing in the ensuing analysis will be assumed to display the same 
noncommutative parameter 6. We shall frequently use the notation: 



a A 6 = a^'Eb = ^ ^ijCLibj = 0162 — 0261. 



Various arguments lead to the following estimate [16]: 



6* < 4 X 10~^°m^ 



(4) 



(5) 



To derive the noncommutative version of a quantum mechanical system one may follow two equivalent 
strategies. In the first case, one considers the usual Hamiltonian where now position dependent functions 
are multiplied by resorting to the Moyal ^^t-product [20]. If one considers e.g. a Hamiltonian of the form 



where q = {qi,q2), = Pi+ P2, then the eigenvalue equation reads: 
where Vq = {d/dqi,d/dq2) and 



(6) 



(7) 



A{q) *e B{q) = exp ( -Vg A V^. ) A{q)B{^] 



(8) 



q'=q 



One may however follow an alternative (but equivalent) route by noting that through a sort of " Seiberg- 
Witten map" [8] 

- 

Ri = qi + Q^^ijPj^ = (9) 

the noncommutative algebra (3) may be transformed into the usual Heisenberg algebra: 



0, 



Ri, Ilj 



(10) 



In eq.(9) and henceforth we shall adopt the Einstein convention. Once the Hamiltonian (6) has been 
expressed in terms of the canonical variables {R, H), the usual quantization procedure follows. In particular, 
the space of states of noncommutative quantum mechanics remains 
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The aim of this research project is to study the emergence of ordinary quantum mechanics in the 
context of noncommutative quantum mechanics. Wc shall follow [33] and call this the noncommutative- 
commutative (NC-C) transition. This question is as pertinent as the question of how the classical world 
emerges in quantum theory [34]. It is well known that, in this case, there is no simple answer. Depending on 
the system [33] the criteria for classicality may vary. One may consider the formal limits ^ — ^ 0, n — > (n is 
some quantum number), N ^ oo (N is the number of particles), or, alternatively, various approximations 
such as Ehrenfest's theorem, WKB approximations, high-temperature expansions or coherent states. In 
many cases one even resorts to combinations of various of these criteria. 

Our strategy to induce a NC-C transition consists of coupling our noncommutative system to an 
external bath and thus treat it as an open system. This is well known in the context of environment- 
induced decoherence where a Brownian particle interacts with a heat bath at thermal equilibrium [34] . To 
keep our calculations as simple as possible we shall assume as in [35] that the reservoir is constituted of 
identical bosonic particles interacting through a potential which has an absolute minimum and that there 
are no appreciable deviations from the equilibrium. Moreover we require that the coupling be weak so 
that only the linear response will be considered. This is commonly known as the harmonic approximation. 
The difference here is that all the particles involved live on the noncommutative plane with identical 
noncommutative parameter 6. In the present paper we shall obtain an exact master equation (to all orders 
in 6) for the reduced Wigner function of the Brownian particle which is the noncommutative extension of 
the Hu-Paz-Zhang equation [36] , [37] . In our derivation we shall follow closely the method of Halliwell and 
Yu [38] . For this purpose a Weyl- Wigner formulation for noncommutative quantum mechanics is required 
[39]- [45]. Some work in this direction has already been developed [8], [33], [46]. For our purposes we just 
need to recapitulate a few basic facts. In addition to the Moyal ^g-product defined in (8) which involves 
only the position variables we shall also consider the usual Moyal ★/j-product: 



(11) 



Aiz) *n Biz) = exp (^^V, ■ JV,'^ Aiz)Biz') 
where = {d/dq, d/dp), z = {q,p), etc and J is the 2d x 2d symplectic matrix: 

(12) 




and d is the number of degrees of freedom (in this paper d = 2). However, the basic algebraic structures 
for noncommutative quantum mechanics in phase space are the extended ★-product (d = 2) : 



Aiz)* Biz) = Aiz) i<n *eBiz) = exp (^y • JV^' + y Vg A V^/) Aiz)Biz') 
and the extended Moyal bracket: 



(13) 



[Aiz), Biz)] = 1 iAiz) * Biz) - Biz) * Aiz)) . (14) 

These expressions were derived in refs.[33, 46, 47] using various methods. The dynamics of the noncom- 
mutative Wigner function is dictated by the extended Moyal equation: 

dW 



dt 



■iz,t) = [Hiz),Wiz,t)]. (15) 



In the previous equation, the noncommutative Wigner function VF(z) describes the state of the system. If 
the state is pure, then (for d = 2) [46], [47]: 

Wiq,p) = -4t2 / , dy e~>PiPiq - y) i^e W+v), (16) 
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where the -kQ product, which acts only on q, is given by (8). As usual, mixed states are described by 
convex combinations of states of the form (16). An important thing to remark is that, in general, these 
noncommutative Wigner functions are not the ordinary Wigner functions of quantum mechanics. Indeed, 
if we compute the marginal distribution 

[ dpW{q,p)=2l;{q)ice^, (17) 

we conclude that it may not be everywhere non-negative. And thus, it cannot be interpreted as a true joint 
probability measure for qi,q2- Clearly, this is a consequence of the noncommutativity (3) and the associated 
uncertainty relations. Nevertheless, the existence of functions which are simultaneously noncommutative 
and ordinary Wigner functions is not precluded [48]. In this work, we shall use as a criterion for the 
NC-C transition, the fact that the quasi-distribution describing the state of the Brownian particle is both 
a noncommutative and an ordinary Wigner function for a sufficiently long time (typically the relaxation 
time scale on which the particle reaches equilibrium with its environment). 

The noncommutative Wigner function for a state at equilibrium (/3 = -j^^) with density matrix 



(Xre 



is given by [46] 

W{z) = Ne-f^^, (19) 
where is a normalization constant, z may this time represent a collection of positions and momenta for 

/3TT 

several particles and ^(/5, z) = is the noncommutative exponential solution of the equation (regarded 

as a partial differential equation) : 

= + (20) 

and subject to the boundary condition: 

</,(/3 = 0,z) = 1, Vzer*M. (21) 

This is an extension of the formula obtained for the deformation quantization of ordinary quantum me- 
chanics in [49, 50]. In these references the extended ^-product appearing e.g. in eq.(20) is replaced by the 
Moyal product ★/j (11). 

It is also important to note that the quantity 

p[W] = J dz W^iz), (22) 

associated with the state W is still a measure of the purity of the system. In rcf.[46] the authors proved 
that there is an upper bound for the purity just like in the commutative case. Indeed if the system is in a 
state represented by the noncommutative Wigner function W{q,p) then: 

p[W] = ^2nh)'^ ' W{z) is a pure state 

(23) 

p[W] < ^2^f^)<i ' W{z) is a mixed state 

This is all that will be required for the purposes of this paper in terms of the Weyl- Wigner formulation of 
noncommutative quantum systems. We have addressed this issue in more depth elsewhere [47]. 
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Before we conclude, a few remarks are in order. It may seem unnatural to have spatial noncommuta- 
tivity and not the full phase space noncommutativity with the last commutation relation in (3) replaced 
by [15], [30]: 

[Pi,Pj] ='iV(^ij, «,J = 1,2. (24) 

We shall however consider the simpler version in (3). There are various reasons for this, (i) Most authors 
[20], [31], [22] regard the algebra in (3) as the basis for noncommutative quantum mechanics. This is fact 
a direct prescription stemming from string theory, where the momentum noncommutativity is absent, (ii) 
Moreover, there are some interesting connections with the Landau problem and the quantum Hall effect 
[31], [51], [52]. (iii) For the particular physical situation of the present work, the derivation of the master 
equation is cumbersome and an additional deformation of the momentum commutation relations would 
make this task even harder, while undermining our main objective: that of showing that the dissipative 
interaction of the system with an external environment is a suitable mechanism to induce the appearance 
of ordinary quantum mechanics in the realm of noncommutative quantum mechanics (we shall address 
this issue in section 6). (iv) As we mentioned before, in order to assess whether the NC-C transition has 
occurred, we need to know whether the state of the Brownian particle is described by a quasi-distribution 
which is simultaneously a noncommutative and an ordinary Wigner function. While it is by now firmly 
established what is meant by a noncommutative Wigner function (cf.(16)), when there is only spatial 
noncommutativity, there is no undisputed counterpart when the momenta are also noncommuting. For 
these reasons, we shall stick to the algebra (3). 

Finally, let us summarize our motivation for writing this paper. Our main purpose is to derive a master 
equation governing the behaviour of a Brownian particle in interaction with an external environment in 
the context of a noncommutative extension of quantum mechanics. This equation provides an interesting 
starting point to address various (conceptual and structural) open issues on noncommutative quantum 
mechanics, namely: 

(i) Can we understand under which conditions is a noncommutative system accurately described by the 
rules of ordinary quantum mechanics? This is an important point: if our universe is noncommutative 
and noncommutative quantum field theories are thus more fundamental than ordinary quantum theory, 
then we should be able to explain how ordinary quantum mechanics emerges in this context. The lessons 
learned in the past decade about the transition from ordinary quantum mechanics to classical mechanics 
will hopefully be a useful guide. One approach that has been intensively explored in this context is that 
of dechoerence. Can a similar mechanism of dechoerence explain why our world looks "commutative" just 
as it partially explains why it looks mainly classical? In this paper we prove that this seems to be the 
case. In fact the example studied in section 6 suggests that, in agreement with the claims of [33], the 
noncommutative-commutative transition takes place before the quantum to classical transition. Moreover, 
the time scale of the former transition is (at least) 7 orders of magnitude shorter than that of the latter. 
And both of them are extremely short. This highlights the fact that we see a classical world and with some 
effort a quantum world, but not a "noncommutative" quantum world. 

(ii) Another important aspect concerns the symmetries of the theory. Many theoretical and experimental 
works (see e.g. [53]-[58]) have recently addressed the breakdown of Lorentz invariance. Noncommutative 
quantum field theories are known to break this symmetry [11]. The system considered in this work is 
non-relativistic. However, we may regard it as a simplified toy model for addressing this problem. Indeed a 
two-dimensional harmonic oscillator Hamiltonian commutes with the angular momentum L = qip2 — q2Pi , 
which is the generator of rotations on the plane. Once the oscillator is coupled to an external environment, 
then its angular momentum ceases to be conserved, and the rate of change is then determined by the torque 
exerted by the environment. Nevertheless this effect is controllable, and if the coupling is sufficiently weak, 
then during a finite lapse of time the angular momentum is approximately conserved. If we now add the 
extra spatial noncommutativity (3) , then the rotational invariance is automatically broken even if we switch 
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off the coupling to the environment. In this noncommutative case, there are thus two contributions to the 
breakdown of rotational invariancc: the torque of the external bath, and the noncommutative contribution. 
If the coupling is sufficiently weak, then during a certain lapse of time (say ti), the external torque will 
have only a residual influence on the angular momentum of the Brownian particle. However, as mentioned 
in (i), the time scale to for the noncommutative-commutative transition is extremely short. After this 
instant Iq, we expect to see the noncommutative contribution to the angular momentum balance equation 
gradually fading away. And thus during the time interval [toi^i]) a rough conservation of the angular 
momentum is to be expected. This would be an interesting toy model to explain why Lorentz invariance is 
such a (experimentally) robust symmetry (even if our universe is noncommutative). We hope to consider 
this problem in a future work. 

This paper is organized as follows. In the next section we obtain an expression for the noncommutative 
Wigner distribution of a collection of noninteracting noncommutative harmonic oscillators at thermal 
equilibrium as well as various expectation values which will be useful for the sequel. In section 3 we derive 
the exact noncommutative Hu-Paz-Zhang equation. In section 4, we obtain the equation that governs the 
dynamics of the purity and state sufficient conditions for a strict decrease of this quantity. As particular 
cases we address the weak coupling limit in section 5, and a noncommutative version of the Caldeira- 
Leggett model [35] in section 6. In the latter case, we obtain an estimate for the time scale of the NC-C 
transition. Finally, in section 7 we present our conclusions. The lengthy derivations of the coefficients of 
the master equation and of the equilibrium distribution of the bath are relegated to the Appendices. 



2 Equilibrium distribution and momenta for a heat bath of noncom- 
mutative harmonic oscillators 

The purpose of this section is to state useful results concerning the equilibrium noncommutative Wigner 
distribution for a bath of N noncommutative harmonic oscillators at temperature T = All the 

formulae presented here are derived explicitly in Appendix 1. We assume the oscillators to be noninter- 
acting, and to have masses (m„) and frequencies (a;„). The corresponding variables ^g^^^p^"^^ obey the 
noncommutative algebra: 



The Hamiltonian is: 



Jn) Jm) 

Qi ,Pj 



-s(rt) Am) 
Pi ,Pj 



0. 



N fn(")V 1 9 



n=l 



The following quantities will be useful: 



mr. 



As the oscillators are noninteracting, the noncommutative Wigner function factorizes: 



(25) 



2mn 



, Vtn = WnA/l + (AnM)2. (27) 



(^tanh[M(j)^ + A„) 



tanh 



2 i^n ^nj 



exp 



-an(/3) - c„(/3) (gW)' - 26„(/?)lW 

(28) 
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where: 



^niP) — ^2hMn^i ^^^^ ^{^n " \ 

Cn(/3) = ^{tanh[M(Q^_A„) 



(0„ + A„)] } 



+ tanh 



bniP) 



(Qn+An) f„„X, 



2^ (^n Afi,; 



^^^^f^tanh ^(Qn + A„ 



For the purposes of this work, we shaU need the foUowing expectation values: 



< ql ' >=< p\ >= 



^ 2 4n„ 



, An) 4m) _ n5„,mSijM„n„ 



^ 2 



{coth (J)^ + An)] + coth [M (o„ - An)] } 



(29) 



[{fin + An) coth (J^n + An)] - (J^n - An) COth (J^^ _ An)] } 



>= {(f^n + An)' coth [M (J^^ + An)] + (J^n - An)' COth [f (J^^ _ An)] } 



(30) 



3 The noncommutative Brownian particle 

In this section we shah consider the combined system of a noncommutative quantum Brownian particle 
and the heat bath of the previous section. The Brownian particle is an oscillator of mass M and bare 
frequency J7. Its coordinates and momenta are {q,p). We assume the particle to be linearly coupled to the 
environment so that the total Hamiltonian reads: 



-2 -\ N N 



(31) 



n=l 



n=l 



where (") is as in eq. (26) . As usual [35] , [38] we assume that at time t = the system and environment 
are uncorrelated: 

Wo {q,p; {g^'^^P^")}) = Wi{q,p) ■ {{q^''\p^''^}) , (32) 

where Wq is given by (27-29). The Wigner function of the combined ensemble satisfies the noncommutative 
Moyal equation (15): 

in^ = Hi.W-Wi.H^^ = -^.VgW + Mn^q . VpW + En (-^ • V^(„) + mnu^q^^'^ ■ V^w) W+ 
+ En Cn -Vp + q- Vp(„)) + f En ("^n^^^gH + Cnq) A V^(„) W + ^ (mO^^ + En Cnq^""^) A VqW. 

(33) 

The reduced Wigner function of the Brownian particle is obtained from W by tracing out the environment's 
degrees of freedom: 



Wr{q,p) = J 



W(q,p,{q(^\p(^^}). 



(34) 
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Let us now perform this integration in eq.(33). We shall assume that W and its first derivatives vanish at 
infinity, so that there will be no phase-space surface contributions. The result is: 



+^Mn\ A VgWr - f E„ CnVg A / [Um dq^'^Up^"'^] gWW^ {q,p, . 
So we have to evaluate the following quantities: 



+ 



(35) 



Giq,p) = Y,Cn f 



Following the same argument as in [38] we conclude that: 

G{q,p) = [A{t)q + B(t)p + C(t)V, + D(t)Vp] Wr{q,p), 



(36) 



(37) 



where A, B, C, D are, for the time being, arbitrary time dependent 2x2 matrices. We then get from 
(35): 

^ = -M • "^l^r + Mn'^q ■ VpWr + {VpWr) ■ A{t)q + Vp ■ {B{t)pWr) + Vp ■ {C{t)VgWr) + Vp ■ (D{t)VpWr) + 

+ 1 MO^g A VgWr - f A iA{t)qWr) - f A (B(t)pW,) - f A (C(t)VqVF,) - f A (D(t)VpVF,) . 

(38) 

It is important to remark that, in the commutative limit {6 — > 0) and for quadratic Hamiltonians, there 
is no difference between the Moyal equation (15) and the classical Liouville equation. In that case the 
quantum effects manifest themselves through the bath's equilibrium distribution (28) (with 9 = 0). This is 
where Planck's constant makes its appearance. At the level of the master equation these quantum effects 
will be encoded solely in the matrix coefficients A(t), B(t), C(i), D(i). In contrast with this behaviour, in 
the noncommutative case (^ > 0), the quantum noncommutative effects appear explicitly in the master 
equation (38) even for quadratic Hamiltonians. 

We conclude this section by presenting an expression for the matrix coefficients A(t), B(t), C(t), D(t) 
in terms of physical quantities which characterize the environment. All the relevant equations are derived 
in the lengthy Appendix 2. 

As in the commutative case the environment is characterized by a dissipation {r]ij{t)) and a noise kernel 
(I'ijit)). These read: 

^ r+oo 

[l~^{co)6ij cos{Lot) + I~{uj)eij sin(a;t)] , 



Vijit) 



dt 



CO 



and 



r+co 

^iji^) — I du coth 
Jo 

Here /^(o;) are the spectral densities: 



[I^{u!)dij cos{ujt) + / {co)eij sin(a;i)] 



[6{lo - fi„ - A„) ± 5{uj - 0„ + A„)] , 



(39) 
(40) 

(41) 



where A„,f2„ are as in eq.(27). The matrix coefficients of the master equation, will also depend on the 
elementary functions Uij{s), Vij{s), linearly independent solutions of the homogeneous integrodifferential 
equations [36], [37]: 

tij{s) + n^^ijis) - ^Mn^eiktkjis) + ^ [ dT rjki{s - r) (Sik - ^eik^) S,j(t) = 0, (42) 



M 



h dr 
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with the boundary conditions: 



Uij{s = 0) = Sij, Uij{s = t) = 0, 
Vij{s = 0) = 0, Vij{s = t) = Sij 



(43) 



We can subsequently construct the functions: 

Lik{s,T)L'^j^{T,T), if s > r 
0, otherwise ' 



Lik{s,T)L^j {t,t), if r > s 



0, 



otherwise 



with 



o 

Lij{s, r) = Uij{s) - Vik{s) {v{t))^1 u{T)ij, L^^(r, r) = g^Lijis, r) 



(44) 
(45) 



Functions G\f{s,T) {k = 1,2) are the Green functions of eq.(147) (Appendix 2) with boundary conditions: 



Sij, 



Sij, 



(46) 



Here and henceforth, we shall tacitly assume that the inverse matrices of {vij{s)) and {uij{s)) exist. We 
are now in a position to present the expressions for all the coefficients appearing in the master equation. 
The result is: 

Aij{t) = 2 /* ds riki{t - s) [5,k - ^eik£) [vij{s) - u,„(s) {u{t))-J v,j{t) 

Bij{t) = jj /o ds vik{t - s) (Sii - ^eii£) {uka{s) {u{t))-} + '-ferj [vUs) - Uka{s) {u{t))-J nS)] } 

Ci3{t)= ^J^dX Af^(t,X)i.ik{t-X) + ^J^dsf^dTj^dX Af^{t,X)A\^\s,T)r,iiit-s)urk{r-X) 

D,,it)= nj^dX [f^Gf^(t,Xp,kit-X) + '^J^dsj'^dTj^dX r^^i{t-s)[iGf^{t,X)]A^\s,r)i.rkir-X) 

(47) 

(48) 



where 



^S\s,r) = (Su - ^Men^^ G«(.,r), A; = 1,2. 



Two remarks are now in order. First of all the integrals in the previous expressions should be regarded as 
improper in the sense that e.g /q ds = limg^o+ /o^^ ds. This is also required in the commutative case [38]. 
Secondly, as in the commutative case, all the matrices in the previous equation are in principle completely 
determined once the dissipation and noise kernels rjijit), i'ij{t) are given. In particular, they do not depend 
on the initial conditions of the Brownian particle. 

The structure of equations (47) looks very complicated. However, in most applications, one considers 
simplified situations, such as the weak coupling limit (see section 5) or the Ohmic high temperature 
Caldeira-Leggett model (see section 6). 



4 Purity 

A common issue in the context of dissipative quantum mechanics is that of finding general conditions of 
dissipativity [59], [60], i.e conditions under which the purity of a state is guaranteed to decrease. Using 
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the noncommutative Hu-Paz-Zhang equation derived in section 3 we may now obtain a set of sufficient 
conditions for this "H-theorem" to hold for the reduced Wigner function Wr- After a simple calculation 
we obtain: 

p[Wr] =Tr(B- ^EA^ p[Wr] - J dz V,Wr ■ JV.Wr, (49) 

where z and Vz are as in cq.(ll). The 2x2 matrices A, B, C and D are as in (47), E is the 2x2 matrix 
with entries e^j and J is the 4x4 matrix: 



J 



f(C^E-EC) C^-fED 
D + D^ 

Since p[Wr] > 0, we then have the following set of sufficient conditions for strict "purity" decrease: 



C + f D^E 



(50) 



(i) the 2x2 matrix B — f EA has nonpositive trace; and 

(ii) the 4x4 matrix J7 is positive semidefinite. 

5 The weak coupling limit 

As a particular case we consider a weak coupling between the Brownian particle and the heat bath. We 
will determine the matrices (47) to order O (C^). As r]ij{t), i'ij{t) are of order O (C^), we conclude from 
(47) that we shall need the functions Uij{t), Vij(t) only to zero-th order: 



±ij{s) + J^^Sy (s) - -MQ^eik±kj{s) = 0. 



(51) 

T.ij{s) = mikPkj {s) + nik\kj{s), (52) 
' ^ij(^) = ^ T,a=± i^iji^ - f^-^) COS ((E + aX)s) + aeij{T, - aX) sin ((E + aX)s)] 



The general solution is given by: 



where real constants and 



Pij{s) = ^ Ef7=± [^ij sin + (^A)s) - aeij cos ((S + aX)s)\ 



with 



Notice that: 



A = 



2h ' 



1 + 



n 



(53) 



(54) 



Xij (s = 0) = pij (s = 0) = 5ij, Xij (s = 0) = pij (s = 0) = 0. (55) 
To compute Aij{t), Bij{t) we do not really need the functions Uij{s), %(«), but rather the alternative set 

Rij{s) = Vij{s) - Uik{s) (n(t))fc/ vij{t), Tij{s) = Uik{s) (n(t))^/ , (56) 

which are equally linearly independent solutions of (42) with initial conditions: 

Rij {s = t)= tij {s = t)= dij , Rij {s = t)= Tij {s = t)=0 (57) 

Prom (52,55,57) we obtain: 

Rij{s) = Xij{s - 1), Tij{s) = pij{s - t) (58) 
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We then get from (47): 
< 

In the commutative hmit (A — > 0, E — > il), 



(59) 



Pjii 



0M 

h ""3 



where: 



Xijit) dij cos(rit), pij{t) sin(Ot), riij{t) Sijr]{t), 



Cl 



MiOn Jo 



ds ri{s) sm{ns) 



is the commutative dissipation kernel. We then have: 

Aij{t)^2Sij [ ds 77(5) cos(Os), Bij(t) 
Jo 

as expected [38]. To compute Cij{t), we need (s < r): 

Gif{s,T) = [Uik{s)vkl{T) - Vik{s)ukl{T)] [uir{T)Vrj{T) - Vlr{T)Urj{T)] 

Again, this is a solution of (42) in the variable s with conditions: 



-1 



0, 



ds ^ ' ' 



Jij- 



Therefore the solution is: 



The same holds for G^- (s,t) with s < t. Notice that in eq.(47) the second term in the expression for 



Gf^{s,T)=Pi,{s-T). 



(60) 
(61) 

(62) 

(63) 

(64) 
(65) 



Cij{t) is of order greater than O (C^). We are thus left with: 



Similarly: 



In the commutative limit: 



^'^■^^^ ^Jl Jo^^ ^^'^^^^ ('^^■' ~ ^^^■^^) 
Dij(t) = n J^ds Uik[s) [dji - -^^ji-j^ J Plk{s)- 



where 



i^is) 



Cij{t) ^ —jpr / ds z/(s) sin(ns), Dij{t) ^ HSij / ds i^(s) cos(Os), 
Mu Jo Jo 



f fil3uj\ 

j div coth y 2 J ^^'^^ cos{u)t), = ^2 



CI 



6{UJ - LUn) 



(66) 
(67) 

(68) 
(69) 



are the commutative noise kernel and spectral density, respectively. Again this is in agreement with known 
results. 
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6 The Caldeira-Leggett model 



In this section we consider a simplified situation known as the Caldeira-Leggett model [35]. We shall 
assume a continuum of oscillators with typical mass: 

rUn ~ m, Vn (70) 

such that: 

r+oo 

Y^d^ du PDiL0)C\L0), (71) 

where 



,0; < A 



PDiio)C\u;) = \ (72) 
[0 ,a; > A 

Here A is a high-frequency cutoff, and r/ is a damping constant. This model is Ohmic in the commutative 
limit in the sense that the spectral density (41) /"""(w) is proportional to a;. In fact: 

7+(a;) = ^, 7-(a;) = (^ = 0). (73) 

TT 

We shall now make several approximations. We shall assume a high-temperature regime leading to a rapid 
localization, so that the localization timescale is much shorter than the relaxation timescale on which 
the system reaches equilibrium with its environment. Moreover we shall only consider noncommutative 

corrections to first order in 6. This means that our description of the system will apply to short times and 
that is small when compared with the typical coherence length. Altogether: 

Finally, we consider the typical mass of the bath to be much smaller than that of the Brownian particle: 

Under these assumptions and after some calculation, we obtain from (39,40): 

%(i)~ry(5., + nfe.,|) 5'it) 

(76) 

^^J{t) - {Sij + S{t) 
where as usual, we used the representation of the delta function: 

S{t)= hm (77) 

Here and henceforth, the symbol ~ denotes the approximation (74,75). 

To compute the coefficients A, B,C,D of the master equation, we need the functions Rij{s),Tij{s) as 
they were defined in (56). In fact, we shall only need to know their boundary conditions (57), and their 
differential equation (42) to order zero in 9. Prom (47,57,76), we obtain after some integrations by parts: 

Mt) - -^euRi^\t) - 2r^5^m - ^S'^j, (78) 
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where rIj\s) is the zeroth-order (in 9) approximation of Rij{s). To this order, we have from eq.(42): 

ii^Ht) = -J^?en<(0 - = -^reJij. (79) 

Here ^ren is the renormahzed frequency [35], [38]: 

^ren = ' ^S{0). (80) 

Altogether: 

Aij{t) ~ -27j6^j6i0) + ^^^72^„e,,-. (81) 
Notice that, as in the commutative case, Aij{t) is time-independent. In a similar fashion, we obtain: 

+ ^^^^ 

where we used 

T§\t) = -^S.j. (83) 

Similarly, to compute Cij{t), we perform a few integrations by parts and use the boundary conditions (46) 
for the Green functions G^j^(s,r), g\^\s,t). The result reads: 

Cij{t) - 2r]KBT^eij. (84) 



Notice that this term is absent in the commutative 
Finally, to compute Dij{t), we need '^Gj^\t,T] 
eq.(147) (Appendix 2): 



imit. 

to zero-th order in 9. This can be evaluated from 

T=t 



.(1) 

We then obtain 



T=t 



Dij{t) ~ 27]kBT6ij + 2r]kBT-ae^j, a = 2r] - y<5(0). (86) 

As in the commutative limit, in the Caldeira-Lcggctt model all the coefficients in the master equation are 
constant. Upon substitution of these coeficients in the master equation (38), we obtain: 



_2er^!^^ A V,Wr + ^Vp A {pWr) - ^^V, A VpWr+ (87) 

+ '-^q A VgWr - ^V, A (pWr) 
For this model, we obtain from (87) or (49), the following equation for the purity of the system: 

P [Wr] = [Wr] - 4nkBT j dqj dp | VpW,(g,p)|2 . (88) 

We conclude that the purity obeys the same dynamical equation as in the commutative limit. To observe 
differences in this quantity, we should consider higher order corrections in our approximation. 
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We now solve the master equation (87), but not in its full generality. We shall set f^r-en = and, in the 
spirit of (74), we shall take the limit T — +oo, ?7 — 0, while keeping the product: 

r = 2r]kBT (89) 

constant. Altogether, we obtain: 

^ ~ • ^^^'^ + ^^i^^ - ^ ^^^'^ ^^^^ 

Strictly speaking this is not a Fokker-Planck equation, as the diffusion matrix is not positive defined. 
Similarly, in eq.(88) we obtain: 

p [Wr] = -2T Jdqjdp \VpWr{q,p)f . (91) 

The purity is thus strictly decreasing. 

Using standard methods [61], [34], we obtain the solution: 

Wr{q,p, t)^ j dq' J dp' Gt (q-^P- <i',P- p') Wo{q',p'), (92) 
where Wo{q,p) is the initial noncommutative Wigner quasi-distribution: 

Wr{q,p,t = 0)=Woiq,p), (93) 

and Gt{q,p) is the Gaussian kernel: 

Gt{z) ~ 3 (^^y exp (-z^Mtz) , z'^ = {q,p). (94) 

The matrix Mj reads: 

3/lf2 / ^2x2 igl2x2 + fE 

M* = ^ . I (95) 



t T e-cr _*?_T 



So Wr{q,p,t) is the convolution (modulo a symplcctic transformation) of a noncommutative Wigner func- 
tion with a Gaussian. Once the Gaussian satisfies the Heisenberg uncertainty relations 

Q 

AqiAq2 ~ - (96) 

this convolution will become a Wigner quasi-distribution [48]. As mentioned in the introduction, we assume 
that for the NC-C transition to occur, then a necessary (possibly not sufficient) condition should be the fact 
that the noncommutative Wigner function describing the state of the system become equally an ordinary 
Wigner function. This seems to be a reasonable criterion in analogy with the transition from ordinary 
quantum mechanics to classical mechanics. Indeed, in this transition the density matrix of the Brownian 
particle becomes approximately diagonal in the position basis. A necessary (but not sufficient) condition for 
this is that the corresponding Wigner function be positive. Notice also that if the noncommutative Wigner 
function becomes equally an ordinary Wigner function, then in particular, the marginal distribution (17), 
will be everywhere non-negative, regardless of whether it was positive or not at the initial time. This will 
allow the usual quantum mechanical probabilistic interpretation for the position variables compatible with 
a course graining resolving areas greater than 9. 
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If we compute A^j for the Gaussian Gt{q-,p)-, we obtain: 



A..-V3M^. (97) 



From (96,97), we thus get: 



(98) 



To this order in 9 the noncommutative Wigner function becomes a Wigner function at te irrespective of 
its initial distribution Wo(q',p). Strictly speaking, this may happen even before Iq. But what our analysis 
shows is that from onwards, it will certainly be a Wigner function. This situation is analogous to that 
described in [62], [61]. 

It is instructive to compare this time scale with the typical time scale of decoherence, where the 
transition from ordinary quantum mechanics to classical mechanics takes place: 



tu ^ r-^. (99) 



Prom (98,99), we obtain: 



Using the estimate (5), we get 



3 MOtl, , , 



< 10-2.^/Mt^, (101) 



if M is expressed in kg and t^, to in seconds. As an example consider an electron in a medium constituted 
of air molecules at T = 300K and pressure of 1 atm. For this situation we have [34] : 

r ~ \^-^'^n?kg^s-^. (102) 

We thus obtain: 

Id ~ 10~^^s te ~ 10~^^s, (103) 

were we used the estimate (5) for d. And thus the time scale for the NC-C transition is much shorter than 
tj). And one should keep in mind that (5) is an upper bound. Therefore, tg should be even smaller than 
our estimate (103). Notice that this result gives countenance to the claim in [33], that the NC-C transition 
takes place before the transition from ordinary quantum mechanics to classical mechanics. 



7 Conclusions and outlook 

Let us restate our results. We computed the equilibrium Wigner distribution for the reservoir of noncom- 
mutative harmonic oscillators as well as the first and second order moments of this distribution. Assuming 
uncorrelated but otherwise arbitrary initial distributions for the Brownian particle and the bath, we derived 
the exact noncommutative extension of the Hu-Paz-Zhang master equation to all orders in Q. By construc- 
tion our equation is Markovian and of the Lindblad form, which ensures the positivity of the evolution. 
Moreover, and contrary to what happens in the commutative limit, the quantum effects appear explicitly 
in the master equation (38) and not just at the level of the matrix coefficients A(t), B(t), C(t), D(t). We 
therefore expect these noncommutative corrections to lead to qualitatively different results. We also stated 
sufficient conditions for strict decrease of purity of the Brownian particle. Finally, we considered the par- 
ticular cases of the weak coupling limit and the Caldeira-Leggett model. In the latter case, we solved the 
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master equation and obtained an estimate for the time scale of the NC-C transition. For an electron at 
room temperature and pressure it is at least seven orders of magnitude shorter than the decoherence time 
scale. Just as the decoherence mechanism explains why universe looks mainly classical, it also seems to 
explain why it is mainly commutative. Our result supports the claim made in [33] that the NC-C transition 
takes place before the quantum-classical transition. 

In a future work, we expect to further explore this master equation [48]. In particular, we intend to 
study a coloured environment, entanglement, and investigate the appearance of certain symmetries once 
the NC-C transition occurs. We are obviously thinking of toy models for Lorentz invariance. As there 
seem to be different time scales for the various transitions, we could therefore consider a situation where 
various systems (with different decohering time scales coexist). We may thus envisage the possibility 
of extending hybrid ensembles of coupled classical and quantum subsystems [63]- [66] to hybrid classical- 
quantum commutative-quantum noncommutative systems. 



Appendix 1 

Our aim is to derive the various formulae presented in section 2. Since the oscillators are noninteracting, 
we may start by evaluating the distribution for a single particle of mass m and frequency uj. Its position 
and momentum variables are q = (^1,^2), P = {pi,P2), respectively, and we assume that they obey the 
noncommutative algebra (3). The Hamiltonian reads: 



(104) 



In terms of R and ft (cf.(9)), the Hamiltonian may be expressed as: 



H{q,p)=n{R,fl)=^ 



1 + 



/ muj6 



V 2n 



+ -mu^R^ —jO, 

2 2n ' 



(105) 



where C = R All is the angular momentum. The equilibrium density matrix at temperature T is then: 



p = Mexp \^-f3H {q, j5)] = M exp [-/^W (i? , ft) 



(106) 



where is a normalization constant. To obtain the corresponding Wigner function in noncommutative 
phase-space we may resort to the generalized Weyl- Wigner map VF^^^'n) introduced in [39]. Indeed, it is 
easy to prove that the application of this map induces in the phase space T*M the extended T*r-product 
(★ = -kfi-kg) mentioned in the introduction [47]. This map consists of applying the Weyl- map based on the 
Heisenberg algebra (i?, 11) (cf(lO)) and. subsequently, performing the phase-space version of the operator 
diffeomorphism (9). We are thus looking for the noncommutative exponential 



.(/3,i?,n)^exp, /-A 



1 + 



~2n 



- -moj^R^ + -—£ 

2 2h 



(107) 



where -kji is the usual Moyal product based on the Heisenberg variables {R, H) (cf.(ll)). Let us rewrite the 
previous equation as: 



. {P, i?, H) = V {ki, k2, R, H) = exp. 



iki (^ + ^^^'^') +ik2nC 



(108) 
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where 



^ - l+(A/a;)^ : 



ki = i/3, k2 



2% 



(109) 



The exponential (108) is well known. We have computed it elsewhere [40] (see also [41]) in the context of 
the two-dimensional commutative harmonic oscillator. The solution is: 



■0 (fci, /c2, -R, n) = cos ^ 



Ml 
2 



(fel + k2 



cos 



Ml 
2 



(fcl - k2) 



X exp { 4 tan [f (fci + ^2)] + '^Mn^R^ + ^JO.) + ^ tan [f {ki - fe)] + ^M^^B? -^c)]. 

(110) 

To finish our calculation we need to go back to the coordinates (g, p) using (9) . A simple calculation leads 
to: 



nf3 



4>iP,Q,P) = cosh ^ 
where this time L = q Ap and 



cosh 



hp 



X exp 



-a(/3)/ - c(/3)g2 _ 26(/3)L 



(111) 



tanh 



(n-x) 



(0 + A) 



c(/3) = ^ {tanh [^(Q _ A)] + tanh [^(O + A)] } 



(112) 



m = ^ tanh 



hi3 



{n-x) 



i^tanh 



mn + x) 



To check our result, we remark the following. The exponential (111) is formally a solution of (20,21) with 
T*M ~ K^. We have checked explicitly that our expression (111) is indeed a solution. Finally, from 
(19,111), we conclude that the equilibrium Wigner distribution is: 



W{q,p) = 7V(/?)exp [-a{(3)p'' - c{l3)q^ - 2h{(3)L\ 



(113) 



where iV(/3) is the normalization constant. To derive N[f3) let us first define the following generating 
function: 



Z{i,r]) = j dQ j dp W{q,p)e:K.p{^-q + r]-p) 



TT^N 

-62 



ac 



exp 



4(62 - ac) 



The normalization is then given by the condition ^(0, 0) = 1. We thus get: 



N{(3) 



ac — 62 



tanh 



nf5 



(n + x) 



tanh 



hp 



i^-x) 



(114) 



(115) 



We finally have: 
W{q,p) = 



1 



(vr/i)' 



tanh 



hp 



(n + x) 



tanh 



np 



{n-x) 



exp 



-a(/?)p2 - c(/3)g2 _ 26(/3)l] . (116) 



Notice that in the commutative limit ^ ^ 0, the following holds: 

( n^Lo, A ^ 0, b{P) 



a{P) 



^tanh(Y), c(/3)^lf^tanh(M!) 



(117) 
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Consequently: 



Wiq,p) 



tanh ^^^^ 



exp 



2 

- — tanh 



hujP\ I 2 2 

— 7:. ^ -ma; g 

2 \2m 2 



(118) 



which is the correct result for a two-dimensional commutative harmonic oscillator. 



Let us now consider N noninteracting harmonic oscillators with masses m„ and frequencies cOn- The corre- 
sponding variables {^q^'^\p^'^^^ obey the noncommutative algebra (25). Our previous analysis reveals that 
the equilibrium noncommutative Wigner functions is given by eqs. (27-29). The corresponding generating 
function is: 



n=l 



4 (^'n - anCn) 



Let us now compute some expectation values using this generating function: 



< 4"" > 



9 



})=0 



(119) 
0. (120) 



(^(") ,,,(") )=o 



Similarly: 



< >= 



We also have 



({^(n),^(n)|) 



6 

On,m ( 



! d2;b ( tin) (n) 



({^("),?j(")})=0 

■^^^ {(^^n + A„) COth [M (i^^ + A„)] - - An) COth (0„ - A„)] } 



and 



4n„ 



< >= 5„ 9 9 



= {coth (ji„ + A„; 



-I- coth 



(^(") )=o 
^ (^n ~ A„)j I 



(121) 

(^(»),T)("))=0 

(122) 
(123) 



Finally: 



< 



.(n) 4m),^{m) ^(n) 



( c{n) „(n 

n,m ^^(„) (n) > '/ 



(5(1) )=0 



(124) 



hSn 



^ {(J2„ + A„)' coth [M (j^„ + A„)] + {Un - Xnf coth [M (J^„ _ A„)] } 
In the commutative limit we recover the expected results [38]: 

coth (^) 



<. Pi Pj ^ 2 



< Z^' ^' > — > 



(125) 



>- 



M coth (^) 
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Appendix 2 

Our task is now to determine the coefficients A, B, C, D. From eq.(33) we have: 



dt ^ "ii ^ — M ^ ^^J ^ 'iJ ^ ^h^V L^n^n ^ 



(126) 



Alternatively, we may perform the calculation using eq.(38) instead: 

i<Qi> = T + i^n^^ij < Qj > +i^ij^jk <Qk> +ieijBjk <Pk> 

d ^ Piqj+qjPi ^ _ 1 ^ ^ ^ ^ )i/fo2 ^ gjgj+gjgj ^ /i , ^ Wk+Mi ^ r , ^ 4jPk+PkQj ^ , ^ 

(127) 

Equating (126) and (127), we obtain the following relations: 

En Cn < ql"^ >= Aij < qj > +Bij < pj > 

(128) 

Mjk = ^ejiNik = Mn6jk 



where: 



[ Njk = EnCn< qfpk > -Aji < Mi+Mk > _B.i < p^pi > +Djk = %fe + Djk 



(129) 



Now let us go back to eq.(38) and consider all the terms that include contributions from the matrices C 
and D: 



dqjdqu 

{m,j - ^ejkNki) - Aij + ^eJknk^] ^ + {N,j - Q.j) £^ + f 6,,- (M,, - A.^) 0^ = (130) 

~\ ^^^j ^ n^Jk'^^kiJ dp.Qq. ''''ijdpidpj h^V^^i-kQq.dq^^, 

where we used (129) and the fact that M and N are symmetric and antisymmetric, respectively. Since the 
matrices M, N do not contribute to our equation it is perfectly consistent to set them to zero. Altogether 
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we have: 



E„ Cn < qt^ >= Aij < qj > +Bij < pj > 



EnCn<qPqk>=Aji< 



qkqi+qiqk 



> +Bji < 



qkpi+piqk 



> -C 



jk 



(131) 



The coefHcients A, B, C, D may now be determined from the previous equation once we have computed 

Jn) ^ 

the expectation vahics < qi >, < q^ qj >, etc of the Heisenberg picture operators by solving the equations 
of motion. Since the Hamiltonian is quadratic we have the following linear solutions: 



m = a,mm + m)pm + (ai;Hmt\^) + ^Sf ^(o)) 



(132) 



for some time dependent coefficients oiij \ \ 6tc. Notice that we have chosen to express q^^"^ in terms 
of qj{t) and Pj{t), instead of qj{0) and Pj{0). This is for later convenience [38]. If we substitute (132) in 
(131) by keeping (30) in mind, we obtain: 



(n,m) 



(rn) , 
«J7 < 



q}-\0) 



(rn) ^ C'^(")'^i""(0)+'jl""Wi4""(0) 



> 



+ 



+b. 



(n,m) 



(m) ^ g(-)(0)p('")(0)+p(-)(0)g(-)(0) (^) 



Pi (0) 



> 



(133) 



Dijit) 



-MEn,n.C^n<| air ^ 



.(m) ^ 



+ 



+6, 



(n,m) 
ik 



dm) . Pi""(0)j,"">(0)+9r' (0)^^(0) ^ , ;(m) 



pt\o) 



> 



Prom (132,133) the problem is solved once we have the solutions of the equations of motion for the 
Heisenberg picture operators. For the sake of simplicity we shall omit the hats over the operators. There 
is no risk of confusion here as all variables appearing in the ensuing analysis are operators. The remaining 
calculations are lengthy and involved. To keep the task as tractable as possible we have chosen to follow 
ref.[38] step-by-step and to adopt their notation. 

To solve the equations of motion it turns out to be more convenient to consider the following set of 
variables: 

/ Qi = qi + ^^ijPj 



n-Pt [ P> -Pi ' 



The rationale for this choice resides in the fact that 



Pi = MQi, p(")=m„QS"^ 



We then get the following set of equations of motion: 



Qi + ^^Qi + CnQ\ ^ — ^M^^CijQj — f E Cn^^ijQ^ 



(n) 



(134) 



(135) 



(136) 



21 



We start by solving the latter equation. The homogeneous equation reads: 



Upon diagonalization we obtain the following solution: 



or + ^iQr - ^rUnule^.Q^;^ = (137) 



QT'd) = >^'mr(S^ + f^'(t)^^^ (138) 



subject to the initial conditions: 

QH(t = 0) = )(0), = 0) = (139) 

The matrices A'^"^ p^^^ are given by: 

^ifi't) = ^ E(T=± (^n - ^^n) COS [(^^n + crA„) t] + creij (On - o-A„) sin [(r2„ + a\n) t]} 



(140) 

'Ij ^ (^) = 2n;r E(t=± {<^y sin [(f^n + crA„) i] - a€ij cos [(ri„ + aXn) t]} 
where A„ are given by (27). It is useful to remark that: 



The matrix A*^"^ is a solution of the differential equation: 

Ajf +^nAlf -2An6,4? = 0. (142) 

Because of relation (141), p^") is a solution of the same equation. Finally, they obey the following initial 
conditions (of. (139)): 

r ASf(i = 0)=<5,, p^\t = 0) = 

(143) 

. A(;)(i = o) = o pS;)(i = o) = <5,, 

A particular solution of the second inhomogeneous equation in (136) is given by: 

"-^/o'"' .g'(*-»)(fa-^^..4)0.W, (144) 
which satisfies the initial conditions: 

g\''\t = 0)=g\''\t = 0) = 0. (145) 
Altogether, the complete solution of the second equation in (136) subject to the conditions (139) is: 



Ql"^(i) = A&)(t)gf ^(0) + p^\t)^i^ + gt\t). (146) 

TTlr] 



If we substitute this expression in the first of eqs.(136) we obtain after some algebra the noncommutative 
Langevin equation: 
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The inhomogeneous term in (147) is: 

/i(t) = -Y.cJ W(t) + 2\.e«pg'(t)l gf (0) + 



1 + ^) "."?(*) ->«4?w' 



(148) 

Equation (147) describes the motion of a noncommutative Brownian particle in a dissipative medium with 
"random" force fi{t) satisfying (cf.(30,40)): 



< m >= 0, < {m, fj{t')} >= huijit - 1'), 



(149) 



where {A, B} = (AB + BA) /2 is the anticommutator. Notice that in the commutative limit 9^0, 
eqs.(41,140,148) reduce to: 



(150) 



/i (0 ^ — Ylin 



rrin OJn 



which are the expected results [38]. Our aim is to solve eq.(147). As in [38], we solve it for two different 
sets of initial conditions^: 

Qi{s = 0) = Qi{0), Qi{s = 0) = (151) 



and 



Qi{s = t) = Qi{t), Qi{s = t) 



M ' 
M ' 



(152) 



where t > is any given time. The solution of the homogeneous equation (/j(t) = 0) (eq.(147)) is given 
by: 

^i(O) 



Ui{s) = Uij{s) - Vik{s) (^(0))^^/ uij{0) Qj{0) + Vik{s) (^(0))^./ 



M 



(153) 



The particular solution to eq.(147) with initial conditions Lbi(s = 0) = u;j(s = 0) = 0, can be formally 
written as: 

^'^'^ = i^r^^^ Gg)(s,r)/,(r), (154) 
where G^^\s,t) is the Green function (44). The solution of (147) with the conditions (151) is then: 



Qi{s) = Ui{s) + Oi{s) = [uij{s) - Vikis) {v{0))m + ^iki^) (^(0))^/ ^ 



-1 PjiO) 



- En ^ /o dr G)^' {s, r) I Aj^^ (r) + IXr^e^ip^^^' (r) I (0) 



.(1) 



An) 



(155) 



En W /o dr {s, r) [(l + f ) (r) - ^e,,A£^ (r) 



^As in [38] s is henceforth the variable and t a specific fixed time. This is related with the fact that in eq. (132) we 
expressed ql"\t) in terms of qi{t),pi{t) and not qi{0),pi(0). 
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Likewise, from (46,135,155) we get: 

Pi{s) = M [iiijis) - Vikis) (vmkl ^iji^)] QM + ^ik{s) {vmk} Pji^) 
- En Cn /o dr \iG\f {s, r)l f A^^ (r) + 2A„e,,pL") (r)] Q^^ (0) 



(156) 



EnCnl^dr [iG^\s,r)] [(l + f ) (r) - ^e,A[?(r) 



Next, we look for the solutions of the homogeneous eq.(147) {fi{t) = 0) with initial conditions (152). The 
result is: 



1 1 1 P'i^) 

Zi{s)= Vij{s) - Uik{s) {u{t))^l Vlj{t) Qj{t) +Uik{s){u{t))^^ ^ 



The solution of the inhomogeneous eq.(147) with the conditions Zi{s = t) = Zi{s = t) = 0, is given by: 

1 r . ^(2) 



(157) 



Zi{s) 



M jt 



dr G\;\s,T)fjiT), is<t) 



(158) 



where G^'^\s) is the Green function (44). Altogether, the solution of (147) with the conditions (152) is 
then: 



Qi{s) = Zi{s) + Zi{s) = Vij{s) - Uik{s) (?i(t))fe/ vij{t) Qj{t) + Uik{s) {u{t)),^j ^+ 
+ En % /.* dr G\f {s, r) [\%^ (r) + 2A„6,,pS,") (r)] Q^"^ (0)+ 

+ E. W fsdr G^is,r) [(l + f ) ,g)(.) - ^.,aS?(.)] 
From (46,135,159) we get: 

Piis) = M [vij{s) - iiikis) {u{t))^l vij{t)] Qjit) + iiikis) («(i))fe/ Pjit)+ 

+ En Cn /i dr [£G\f {s, r)] [a^^ (t) + 2Ane,,p;,") (r)] Q^'^ (0)+ 
+ EnCnl!dr [iG^\s,r)] [(l + f ) pg)(r) - ^6,,Al?(r) 



(159) 



(160) 



Substitution of (159) in (144,146) yields: 



Qf^w=Ar(i)Qr(o)+prw 



/O - S) {S,, - ^e,u£) { [v,l{s) - U,a{s) {u{t))ab Mt)] Ql{t) + U,a{s) m)al ^ 



•lit) I 
M "I" 



+ En. ^ 11 dr G^: is, r) [\'-> (r) + 2\^earpT {r)\ QT" (0)+ 
+ E. )(.,.) [(1 + ) ,M(r) - ^e„.AM(.)] 5^} . 



(161) 
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And similarly: 



-Cn Jo* ds pg)(i - S) {Si, - ^e,^£) { [v,l{s) - U^a{s) {U{t))ab ^^K*)] Ql{t) + «,a(^) 



+ J:m^Isdr Gyj{s,r) A^^ (r) + 2A^e„,p^r W ^ 



Jm) 



p/"'(o) 



(162) 



+ E„ % /i Gfj {s,r)[{l + f^) plr^ (r) - ^e^.A^^^ (r) 

Using these expressions as well as the relations (133,134), we obtain after some manipulations: 

= /o ds [(l + -s)- ^6,.Ay (i - .)] (5,, - ^e.,i) [vi,{s) - u^s) {u{t))a,' v,jit) 



X {uka{s) {u{t))J + ^€rj [vkr{s) - Uka{s) {u{t))J Vbr{t)] } 



(163) 



And also: 

(n,m)/,N _ ri,(n). 



a)?"" (t) = \>^' (t) + 2Xne,kf^ {t)\ 6n,m -^I^dsf^dr [A^^ (r) + 2A^6,.p:;;'^ (r 
xAS?(-,r) [(l + f ) ^(t -s)- ^e,A(?(t - .) 



(m) , 



(n) / 



(164) 



1±Mlodsndr pi7^(r)A|f (.,r)x 



Moreover: 



[Ag)(.) + 2A„6fe,p!;)(.)lA«(t,.) 



Finally, if we substitute eqs. (163-165) in (133), we recover (47,48). 
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